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Abstract 



In this article, we study the decays ij}' — > J /^tt'^'k^ and 77^ — > rjc'n^'!^^ by taking 
into account the chiral symmetry breaking effects, the final-state interactions and the 
^ ' heavy quark symmetry. We can confront the predictions of the ry^ — >■ rjc'K^'K~ decay 

width and differential decay width with the experimental data in the future, and obtain 
powerful constraints on the chiral breaking effects and the final-state interactions, and 
test the heavy quark symmetry. 
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^ ! 1 Introduction 



Hadronic transitions among the charmonium and bottomonium states V'(mS) tp{r]S>)'K~^'K~ 
^ \ and T(mS) T(nS)7r^7r^ are of particular interesting for studying the dynamics of both 

the heavy quarkonia and the light mesons. Such processes are usually calculated with the 
. multipole expansion in QCD, where the heavy quarkonia are considered as the compact 

\ and nonrelativistic objects and emit soft gluons which hadronize into the light meson or 

light meson pair pj. The amplitudes can be f actor ized into the heavy quarkonium part 
. and the light meson part. The former part depends on the dynamics of the quarkonium 

and should preserve the heavy quark symmetry, while the latter part depends on the chiral 
dynamics and should obey the chiral symmetry [21 [3l H] . 

In Ref . [5] , Manuel and Urech construct an effective Lagrangian for the hadronic decays 
^ \ if:' — 7> J/V'vr^TT" and T' — )• T7r"''7r~ based on the heavy quark symmetry and the chiral sym- 

metry, and obtain reasonable values for the coupling constants by fitting to the invariant 
7r"'"7r~ mass distributions. In Ref. [6], Yan, Wei and Zhuang observe that there are D- 
wave contributions besides the 5- wave contributions, although the D-wme contributions 
are very small. The D-wave contributions were firstly predicted by the Novikov-Shifman 
model based on the multipole expansion in QCD combined with the chiral symmetry, 
current algebra and partially conserved axial- vector current [7]. In Refs.[8l|9], Guo et al 
observe that the final-state interactions play an important role. In the case of the transi- 
tions T(3S) T(lS)7r7r, there is a double peak in the tttt invariant mass spectrum [TO] , 
we have to resort to additional assumptions to describe the experimental data, such as the 
relativistic corrections |3], the final-state interactions and the /o(600) resonance [TTt [T2]. 
the exotic Tvr resonances [El dU [H] , the coupled channel effects ^4j , the S — D mixing 
[15] , the field correlators [16], etc. 

The two vr transition r]'^ r/c vr+vr^ has not been observed yet. Recently, the CLEO 
collaboration searched for the decay V' ~^ iv'c ^ sample of 25.9 million tp' events 
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collected with the CLEO-c detector, and observed no evidence for the decays V' — ^ 7^c 
and — >■ r/cTT+Tr", and set the upper limit, 

Br{tp' 777^,) X Br(7/; r/cVr+Tr") < 1.7 x 10~^ , (1) 

at the 90% confidence level [T7j. It is interesting to make predictions for the decay width 
and differential decay width of the process rj'^ — >• ricTT~^iT~ , which may be observed at the 
BESIII and PANDA in the future [151 [T9]. 

In this article, we study the decays ^p' — > J/ip-ir~^7r^ and 77^ — > r/c vr^vr^ with a phe- 
nomenological Lagrangian by taking into account the chiral symmetry breaking effects, 
the final-state interactions and the heavy quark symmetry O [6l [H [9l [20l [21] . 

The article is arranged as follows: we study the 7r"^7r~ transitions of the ip' and ry^ in 
details based on the heavy quark symmetry in Sec. 2; in Sec. 3, we present the numerical 
results and discussions; and Sec. 4 is reserved for our conclusions. 

2 The 7r+7r~ transitions of the i/j' and r/^ with the heavy 
quark symmetry 

The charmonium states can be classified according to the notation n^^^^Lj, where the n 
is the radial quantum number, the L is the orbital angular momentum, the s is the spin, 
and the j is the total angular momentum. They have the parity and charge conjugation 
P = (— 1)''""'"^ and C = (—1)^"'''^, respectively. The states have the same radial quantum 
number n and orbital momentum L = can be expressed by the superfield J \20\ [21] , 

J = -^[4^,r-Vcl,]^, (2) 

where the v'^ denotes the four velocity associated to the superfield. We multiply the 
charmonium fields ip^ and rjc with a factor y^M^ and ^JM^^ respectively, and they have 
dimension of mass |. The superfields have been used to construct the phenomenological 
Lagrangians to study the radiative transitions, pseudoscalar meson transitions and vector 
meson transitions among the heavy quarkonia [201 [2H [22] . 

The 7r"^7r~ transitions between the m and n charmonium states can be described by 
the following phenomenological Lagrangian [5l [21] , 

£ = ^ ^ Tr [J(m)j(n)] |5i(m, n)TV[(9„;7)(a"C/)^] + 52("i, n)Tr[(u • • SC/)^'] 

m,n 

+53(m,n)Tr[M(C/ + C/t - 2)]} , (3) 

where J = 7° J'''7'^, M = Bq diag{mtt, m^, rus} with = 2BQm, mj^ = BQ{m+ms), = 
^Bo(m + 2ms) in the isospin symmetry limit = rud = m, the U is a 3x3 matrix that 
contains the pseudoscalar Goldstone fields, and the gi{m, n) denote the coupling constants. 
We construct the chiral symmetry breaking term |g3(m,n)Tr [j{m)J{n)~\ Tr[M(C/ + — 
2)] consulting Ref.[5]. In the following, we will smear the indexes {m,n) in the coupling 
constants gi{m,n) for simplicity. 
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T^^j/^^^ = —72 [aiPi ■ P2 + 92P1P2 + 53^2] Jm^,Mj/^€* ■ e' , 

In 



The tree diagram amplitudes for the ir'^Tr transitions of the iJj' and r]'^ can be written 

as 

4 

72 

J n 

4 

72 

J n 

where the vr decay constant 7,^ = 92MeV, the pi and p2 are the four-momenta of the vr^ 
and 7r~ respectively, the and P2 'ire the energies of the vr"'" and vr" in the laboratory frame 
respectively, and the e and e' are the polarization vectors of the ^/J' and J/V' respectively. 
The Pi and P2 can be written as functions of the momenta of pions in the center of mass 
frame of the vrvr system: 

= r^ {Ei + m\pi\cose:), 

= -^=L=iEi-mpi\cose:), (5) 

where the /3 is the velocity of the vrvr system in the center of mass frame of the initial 
particle, the pl={E^, p^) and P2={E2, P2) are the four-momenta of the 7r+ and vr" in the 
center of mass frame of the vrvr system respectively. The P1P2 can be written as 



P1P2 



1-/32 



(6) 



where |pi| = jpal = y ^ -ml, the Po(cos6'*) = 1 and P2(cos6'*) = ^{cos'^ 9* - i) are 
the Legendre functions [8l |9] . 

In Refs.[8[i9j, Guo et al retain the coupling constants gi and g2, and observe that the 
5-wave vrvr final-state interactions play an important role and should be properly included. 
The chiral unitary theory is a suitable approach for taking into account the infinite series 
of the re-scattering meson loops [23]. At the lowest order, the isospin 7 = kernel of the 
Bethe-Salpeter equation is 

g _ 

^7r,7r7r('*) ~ 72 ' C^) 

Jit 

and the D-wave final-state interactions cannot be included. 

The scattering amplitudes for the decays — )> J/^tt^tt^ and rj'^ — t- r/c vr"''vr~ can be 
written as 

T = T^ + Tl. G{mt) ■ 2ri=o,(mL) , (8) 

where 2T/^^^ = (vr+vr^ + vr^vr+ + vr'^vr°|T^^'^[vr+vr^), the Tg are the 5-wave components 
of the scattering amplitudes T", and the G{p^) is the two-meson loop propagator. 



(2vr)^ — ml + ie {p — qY — + ie 
1 f~. ^ , rnl , cr -I- 1 ' 



(4vr) 



|5(^) + log^ + alog^| , (9) 
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where 



1 



4m|_ 



and n = nip = 770 MeV. Here we take the dimensional 



regulation to regulate the ultraviolet divergence and introduce the subtraction constant 
a(//) as a free parameter. The full 5-wave tttt vrvr scattering amplitude T^^^^ can be 
taken as the solution of the on-shell Bethe-Salpeter equation 



TTTT, TTTT V "TTTT/ 



l-G(m2jy/=0^(m2, 



(10) 



where we have neglected the contributions from the KK channels considering the values 
- Mj/^ < 2m K and M^/^ - M^^ < 2m k ■ 
The differential decay width of the transition ^' — > J/ipir^n' can be written as 



dT 



ip' —^J/lpTr+TT~ 

d/TYti'TrnT 



(2vr)38M2 



^'^\T\^\Pi\Pj/^pdcos{ 



(11) 



where 



Pj/-<P 



Ml, 



{Mj/^ + m^ 



Ml, 



{Mj/^ 



m-n 



2AL 



(12) 



the denotes the average over the polarization vector of the initial state ip' and the 

sum over the polarization vector of the final state J/ip. The corresponding differential 
decay width of the transition r]'^ — t- r/cvr+7r~ can be obtained with a simple replacement. 



3 Numerical results and discussions 

The coupling constants gi, §2 and and the subtraction constant a{^) can be fitted to the 
experimental data on the transition ■(/'' — > J/V-'^r^vr^ from the BES collaboration [24j. In 
Refs.[5l[6l[8l[9], the coupling constant associates with the small m\ is neglected. In this 
article, we retain the coupling constant 53 and fit the parameters to the experimental data 
in four cases: 2CC, 3CC, 2CC + FSI and 3CC + FSI, respectively, where the 2CC denotes 
the two coupling constants gi and 172, the 3CC denotes the three coupling constants 171, 
g2 and 173, the FSI denotes the final-state interactions. In the chiral limit, the Adler zero 
condition can be satisfied. The numerical results are plotted as the number of events versus 
the TTTT invariant momentum m^^, see Fig.l. From the figure, we can see that retaining 
only the coupling constants gi and 32 can lead to rather satisfactory fitting, by adding the 
coupling constant 33 and final-state interactions, even better fittings can be obtained. 

We normalize the BES data using the width V^, = 286 keV and the branching ratio 
Br('0' — > J/'tpTT~^7r~) = 33.6% [25]. The numerical values of the coupling constants are 
shown in Table 1, where the unit of the coupling constants gi, 52 and gs is GeV~^, the 
subtraction constant a(//) is a dimensionless quantity. From the Table, we can see that 
the parameters from the four cases differ from each other remarkably (or significantly), 
the coupling constant g^ and the final-state interactions maybe play an important role, 
and we should take them into account. 

Using the parameters presented in Table 1, we can obtain the decay width and the 
differential decay width of the transition rj'^ — > r/cTr+vr" in the four cases 2CC, 3CC, 
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2CC 


3CC 


2CC±FSI 


3CC±FSI 


91 


0.0873 ± 0.0008 


0.1086 ± 0.0039 


0.0586 ± 0.0088 


0.0468 ± 0.0014 


92 


-0.0258 ±0.0010 


-0.1814 ±0.0247 


-0.0231 ±0.0031 


0.0033 ± 0.0018 


93 




0.5098 ± 0.0787 




-0.0794 ± 0.0062 


a 






-1.0661 ±0.5755 


-1.8160 ±0.1752 


r 

^ TJcTTTT 


229.7+^:3 


uo.ot^^V 


209.1^t^^i 


240.3^^^:^ 


r 


9 OQ+O.OS 




o 1 0+1.28 


2 cn+O.SS 
^•^'-'-0.35 



Table 1: The parameters fitted to the experimental data on the V' ~^ J/'ijj'K'^TT decay, 
the unit of the r/^ — > r/cvr^vr^ decay width is KeV, and T — — '^c^^c^r^ 



VcTTTT r , 

' '- ip'-tj/ipnir 

2CC ± FSI and 3CC ± FSI, which are shown in Table 1 and Fig. 2, respectively. Prom the 
Fig. 2, we can see that the line-shapes of the differential decay width of the transition rj'^ 
r/cVr+TT^ differ from each other significantly in the four cases, although those parameters 
can all give satisfactory descriptions of the ip' J/^/^vr^vr" differential decay width. 

In Ref . [26j . M. B. Voloshin studies the transitions ip' J/ip7r~^TT~ and r]'^ — )• r]ciT^iT~ 
in the framework of the multipole expansion in QCD using the current algebra and the 
trace anomaly in QCD, and obtain the ratio. 



the lower bound is compatible with the upper bound of the present prediction 2.501q'|5 
in the case of the 3CC±FSI. In the cases of the 3CC and 2CC±FSI, the uncertainties of 
the present predictions are too large. We can confront the present predictions with the 
experimental data at the BESIII and PANDA in the future \18\ [T9]. and obtain powerful 
constraints on the chiral breaking effects and the final-state interactions, and test the 
heavy quark symmetry. 



4 Conclusion 

In this article, we study the decays ip' — )• J/^vr+vr^ and r/^ — )• r/cvr+vr^ by taking into 
account the chiral breaking effects, the final-state interactions, and the heavy quark sym- 
metry. We fit the parameters to the experimental data on the ip' — >■ J/V'tt+tt" from the 
EES collaboration, and then take those values to calculate the decay width and differential 
decay width of the transition ry^ — )■ r]cTT~^iT~ , which can be confronted with the experimen- 
tal data in the future, and put powerful constraints on the chiral breaking effects and the 
final-state interactions, and test the heavy quark symmetry. 
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Figure 1 : The number of events versus the tttt invariant mass distribution 171-,^^^ in the decay 
ip' — 7> J / ipTT~^ Tr~ , the normahzition terms are not shown exphcitly. The experimental data 
is taken from the BES collaboration. 
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Figure 2: The differential decay width of the transition rj'^ — > r/cvr~^7r versus the tttt 
invariant mass distribution m^jr- 
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